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The first N-body simulation of interacting galaxies, even producing spiral arms, was per¬ 
formed by Erik Holmberg in Lund (1941), not with a numerical computer, but by his 
arrangement of movable light-bulbs and photocells to measure the luminosity at each bulb 
and thereby estimate the gravitational force. A decade later, and with the arrival of the 
first programable computers, computations of galactic dynamics were performed, which 
were later transferred into a N-body simulation movie. I present here the background de¬ 
tails for this work with a description of the important elements to note in the movie which 
may be retrieved at http: //ttt. astro. su.se/~po. 

The computer BESK 

In the late 1940ies the Swedish Government set up a State board given the task to con¬ 
struct a modern, new generation computer to meet the country’s computing demands. In 
1953 the computer BESK0 (Binary Electronic Sequence Calculator) was completed and 
ready to start. It then had a Williams Tube based, later changed to ferrite-core, direct 
access memory containing a maximum of 2048 numbers or instructions. An addition or 
subtraction was done in 41 microseconds. The operation times for multiplication and di¬ 
vision were 266 and 512 microseconds respectively. Input and output were performed via 
punch tape. 

The computer was used already from the early stages by Ingrid Torgard in Lund and 
Alexander Ollongren in Leiden (1962) for computations of three-dimensional galactic or¬ 
bits, and by Per Olof Lindblad in Stockholm (1960) for two-dimensional galactic N-body 
simulations. One of these early simulations is shown in the accompanying movie, which 
will illustrate some aspects of galactic dynamics. 


'http://en.Wikipedia.org/wiki/BESK 
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The N-body computing program 


From the very beginning it occurred that an interesting task for this computer would be 
to simulate a galaxy by an N-body problem. The memory size of the computer then puts 
a strict upper limit to the number of bodies. The fourth order Runge-Kutta integration 
method used, requires a fair amount of memory space. As far as only the direct access 
memory was used, and provided the bodies were placed in a two-dimensional and strictly 
bi-symmetric configuration, the maximum number of bodies was restricted to N = 116. 

This restricted number of bodies, confined to a plane, cannot realistically simulate the 
mass distribution in an entire galaxy. Thus, the main part of the gravitational force in 
the plane of the galaxy was represented by a fixed rotationally symmetric central force, 
while the free mass points represented a thin disk. For the fixed central force, the force 
derived by Marten Schmidt (1956) for the Galactic System was chosen, and the mass of 
each point was chosen to be 64 million solar masses. Their total mass of the free bodies 
then would be of the order of 4% of the mass of the galaxy. To avoid violent encounters 
between the mass points, a sphere with ‘soft gravity’ surrounding each mass point had to 
be introduced. The radii of these spheres decide the rate of the ongoing increase of velocity 
dispersion occurring throughout the computation. 

The galactic model 

In the case illustrated in the movie, the mass points were placed in three circular rings with 
radii of 2, 4 and 6 kpc and given initial circular rotational velocities. In order to throw this 
system out of balance, the mass points in the middle ring were shifted slightly towards two 
bisymmetric positions along the ring as indicated in the starting sequence of the movie. 
Then the computation is started, keeping a full integration step of 8 million years, each 
step requiring 7.5 minutes of computing time. The clock in the upper right corner of the 
frame indicates the number of years that have elapsed since the start. 

To note in the movie ( http://ttt.astro.su.se/~po ) 

The immediate reaction on the uneven mass distribution in the middle ring is that the 
points trailing the density maxima are accelerated and move outwards, while the points 
leading the maxima are decelerated and move inwards, giving the ring a slight elongated 
shape. In fact two interfering elongated waves are exited, one moving slower than circular 
motion and hence with density maxima at the ends of the major axis and one moving 
faster than circular motion and hence with density maxima at the minor axis. With these 
two waves superimposed, the result is that the middle ring for a while oscillates between 
elongated shape with even mass distribution, and almost circular shape with bipolar mass 
distribution, i,e, it exerts a bar-like perturbation on its surroundings. The outer ring lies 
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outside the co-rotation resonance of this bar-like perturbation of the middle ring. This 
causes the outer ring to strive for an elongated shape with major axis at right angles to 
that of the bar and with density maxima at the ends of its minor axis. By mutual inter¬ 
actions, the amplitudes increase, and after about 350 million years from the start the two 
rings coalesce, the rings break up and from 400 million years form a pair of leading spiral 
arms. 

These leading spiral arms are not stable. The leading parts of the arms are pulled back 
by the gravitational attraction of the arms, which makes them fall inwards and then shoot 
forwards in what Donald Lynden-Bell has called a ‘donkey behaviour’. The trailing parts 
of the arms, on the other hand, are pulled forward, move outwards and lag behind, again 
a ‘donkey behaviour’. In this way, the middle and the outer arms change place, and after 
550 million years now form a pair of trailing arms. It should be noted that this change 
from leading to trailing arms is not due to differential rotation but to self-gravity within 
the arms that turn them inside out. (This is more clearly seen in the printed version of 
the simulation, where point masses from different rings are marked with different symbols.) 

The trailing arms are subsequently drawn out by differential rotation, until they by 750 
years vanish with increasing velocity dispersion. However, now the inner ring in turn has 
acquired an elongated shape acting as a bar, which again creates a trailing spiral structure 
formed mainly by point masses from the middle ring, 900 million years from the start, and 
after a total computer time of about 15 hours. 

The present movie was prepared during a stay at the Leiden Observatory in the spring of 
1961 for the symposium Interstellar matter in galaxies (ed. L. Woltjer, W.A. Benjamin, 
Inc., p. 222, 1962) held in Princeton in April 1961. Today of course, there have appeared 
innumerable, infinitely more sophisticated, simulations of galaxies. 
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